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Abstract 

We perform a comprehensive study of a number of rare charm decays, incorporat¬ 
ing the first evaluation of the QCD corrections to the short distance contributions, as 
well as examining the long range effects. For processes mediated by the c ^ ui^i~ 
transitions, we show that sensitivity to short distance physics exists in kinematic 
regions away from the vector meson resonances that dominate the total rate. In 
particular, we find that D and D pC^t~ are sensitive to non-universal 

soft-breaking effects in the Minimal Supersymmetric Standard Model with R-parity 
conservation. We separately study the sensitivity of these modes to R-parity violat¬ 
ing effects and derive new bounds on R-parity violating couplings. We also obtain 
predictions for these decays within extensions of the Standard Model, including ex¬ 
tensions of the Higgs, gauge and fermion sectors, as well as models of dynamical 
electroweak symmetry breaking. 


1 Introduction 


The remarkable success of the Standard Model (SM) in describing all experimental infor¬ 
mation currently available suggests that the quest for deviations from it should be directed 
either at higher energy scales or at small effects in low energy observables. To the latter 
group belong the sub-percent level precision measurements of electroweak observables at 
LEP and SLD as well as the Tevatron experiments |l| . Tests of the SM through quantum 
corrections have proved to be a powerful tool for probing the high energy scales possibly 
related to electroweak symmetry breaking and the flavor problem. The absence of flavor 
changing neutral currents (FCNC) at tree level in the SM implies that processes involving 
these currents are a primary test of the quantum structure of the theory. Most of the 
attention on FCNC has been focused on processes involving K and B mesons, such as 
and i?S(s) ~ niixing and also on rare decays involving transitions such as 
s —>• s —>• duu, b —>• S 7 , b —>■ si^i~, etc. 

The analogous FCNC processes in the charm sector have received considerably less 
scrutiny. This is perhaps due to the fact that, on general grounds, the SM expectations 
are very small both for mixing as well as for FCNC decays [^,^,|^. For 

instance, there are no large non-decoupling effects arising from a heavy fermion in the 
leading one-loop contributions. This is in sharp contrast with K and B FCNC processes, 
which are affected by the presence of the top quark in loops. In the SM, D meson FCNC 
transitions involve the rather light down-quark sector which translates into an efficient 
GIM cancellation. In many cases, extensions of the SM may upset this suppression and 
give contributions sometimes orders of magnitude larger than the SM. In this paper we 
wish to investigate this possibility. As a hrst step, and in order to establish the exis¬ 
tence of a clean window for the observation of new physics in a given observable in rare 
charm processes, we must compute the SM contribution to such quantities. This is of 
particular importance in this case due to the presence of potentially large long-distance 
contributions which are non-perturbative in essence and therefore non-calculable by ana¬ 
lytical methods. In general the flavor structure of charm FCNC favors the propagation of 
light-quark-states as intermediate states which, if dominant, obscure the more interesting 
short distance contributions that are the true test of the SM. This is the situation in 
mixing and in the c ^ u'f transition . In the case of mixing, although 

the long distance effects seem to dominate over the SM short distance contributions, it is 
still possible that there is a window of one or two orders of magnitude between these and 
the current experimental limit [§]; the predictions of numerous extensions of the SM lie 
in this window [^]. On the other hand, charm radiative decays are completely dominated 
by non-perturbative physics and do not constitute a suitable test of the short distance 
structure of the SM or its extensions. 

In what follows we investigate the potential of rare charm decays to constrain exten- 
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sions of the SM. With the exception of > 77 , we shall concentrate on the non-radiative 
FCNC transitions such as c —^ c —*• uuu which enter in decays like 

D —>• D —>• XuUu, etc. We extensively consider supersymmetry by studying 

the Minimal Supersymmetric SM (MSSM) as well as supersymmetric scenarios allowing 
R-parity violation. We find that rare charm decays are potentially good tests of the 
MSSM and also serve to constrain R-parity violating couplings in kinematic regions away 
from resonances. In charged dilepton modes, this mostly means at low dilepton mass. In 
general, we hnd that this kinematic region, corresponding to large hadronic recoil, is the 
most sensitive for new physics searches. 


The D —»• decays were studied in Ref. |j^ in the SM without QCD corrections. 

More recently the D —> decays were examined in Ref. in the SM and some of 

its extensions, including the MSSM. We compare these predictions with ours, and hnd 
some discrepancies in the SM calculation of the long distance contributions. We also 
emphasize the importance oi D ^ Vt^^~ in the MSSM due to its enhanced sensitivity 
to the electromagnetic dipole moment operator entering in c —>■ M 7 . 


In the next section we calculate the SM short distance contributions including QCD 
corrections and estimate long distance effects for various decay modes. In Section 3 we 
study possible extensions of the SM that might produce signals which fall below current 
experimental limits but above the SM results of Section 2. We summarize and conclude 
in Section 4. 


As a hnal comment, we note the following convention and notation used throughout 
the paper. Many quantities relating to both SM and also new physics are chiral, involving 
projection operators for left-handed (LH) and right-handed (RH) massless fermions. We 
shall employ the notation 


T . _ 1±75 _ 7^(1 ± 75 ) 

J- L,R — ^ , -L L,R — 2 

for scalar projection operators Tlr and vector projection operators Tljj^. The chiral 
projections of fermion field q are thus expressed as 


( 1 l,r = rL,R O' 


( 2 ) 


2 The Standard Model Contributions 

In this section we study the Standard Model contributions to various charm meson rare 
decays. At the time of this writing, there are no reported events of the type we are 
considering. We group the decay modes by their common short distance structure. In 
each case we address both the perturbative short distance amplitude and the effects of 
the non-perturbative long-range propagation of intermediate hadronic states. Due to 
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the non-perturbative nature of the underlying physics, the long distance effects cannot 
be calculated with controlled uncertainties. Therefore we hnd it prudent to generate 
estimates by using several distinct approaches, such as vector meson dominance (VMD) 
for processes with photon emission and/or calculable unitarity contributions. In this way, 
we hope to obtain a reasonable measure of the uncertainty involved in the calculation, 
and at the same time, obtain bounds on the magnitude of long-distance contributions 
which are not overly model dependent. 

2.1 Meson Lepton-antilepton Transitions D 

As we shall discuss, this mode is likely to be observed at forthcoming B and Charm 
factory/accelerator experiments. We start with the calculation of both short and long 
distance contributions to the inclusive rate. We then compute the rates for various exclu¬ 
sive modes. 


2.1.1 The Short Distance Contribution to D ^ 


The short distance contribution is induced at one loop in the SM. It is convenient to use 
an effective description with the W boson and the b-quark being integrated out as their 
thresholds are reached, respectively, in the renormalization group evolution ||12|| , 


'HeS — 


T~(-eS — 


4Gf 

■ 71 “ 

4Gf 

■ 71 “ 


10 


1 : + Y,Ci(.p)OM 

q=d,s,b 2=3 

10 

c!«'(,.)oS«>(,.) + c<^\p)o!f\p) + Y,c:(p)o',{p) 

q=d,s 2=3 


, mb < n < Mw 
, H<mb, ( 3 ) 


with {Oi} being the complete operator basis, {Gi} the corresponding Wilson coefficients 
and fi the renormalization scale; the primed quantities indicate those where the b-quark 
has been eliminated. Note that we must keep all terms of order 1/M^ above the scale 
/i = mfe in this decay as opposed to radiative decays. In Eq. (H), the Wilson coefficients 
contain the dependence on the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements 
Vgqf. As was pointed out in Ref. [Q, the CKM structure of these transitions is drastically 
different from that of the analogous B meson processes. The operators Oi and O 2 are 
explicitly split into their CKM components 


Oi'^ = {ul7t.ql){qh^cl) , , ( 4 ) 


where q = d,s, b, and a, (3 are contracted color indices. The rest of the operator basis is 
dehned in the standard way. The QCD penguin operators are given by 

O3 = {ui-i^ci)Y,{qiYql). Oa = {ui'^^cl)Y,{fLi^qD , 

<? q 
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05 = , Oq = , (5) 

9 9 

the electromagnetic and chromomagnetic dipole operators are 

Oj = Y^^c{uL(JtiuCR)F>^'' , Os = ■^^rriciuLCTf.iyT'^CRjG'^'' , ( 6 ) 

and finally the four-fermion operators coupling directly to the charged leptons are 

2 2 

Og = ^;^{uL-ff,CL){h'"i) , Oio = Y^(“i7MCL)(^7''75^) • (7) 

The matching conditions at p = Mw for the Wilson coefficients of the operators Oi _6 are 

Ct{Mw) = 0, Cs-e{Mw) = 0, Cl{Mw) =-X, , ( 8 ) 


with Xq = V*qVuq- The corresponding conditions for the coefficients of the operators O 7-10 
are 


Ci{Mw) — 

~2i^sF2{Xs) + XbF2{Xb)} 



Cs{Mw) = 

-^{XsD{xs) + XbD{xb)} , 



cPiMw) = 

i: A. 

i=s,{b) 

- (^Fi{xi) + 2C(xi) 

, , 0(xA 
^ 2.2 J 


C[g{Mw) = 

\ 0{xi) 

AX- 2si ■ 


(9) 


i=s,{b) 


In Eqs. @ we define Xi = the functions Fi{x), F 2 {x) and C{x) are those derived 


in Ref. |T^ and the function D{x) was defined in Ref. 


To compute the c —*• rate at leading order, operators in addition to O 7 , O 9 

and Oio must contribute. Even in the absence of the strong interactions, the insertion of 
the operators in a loop would give a contribution sometimes referred to as leading 
order mixing of C 2 with Cg. When the strong interactions are included, further mixing 
of the four-quark operators with O 7 _io occurs. The effect of these QCD corrections in 
the renormalization group (RG) running from M\y down to /i = rric is of particular 
importance in the coefficient determining the c ^ u'y amplitude. As was shown 

in Ref. P], the QCD-induced mixing with dominates CY^{mc)- The fact that the main 
contribution to the c ^ uj amplitude comes from the insertion of four-quark operators 
inducing light-quark loops signals the presence of large long distance effects. This was 
confirmed in Ref. ^ where these non-perturbative contributions were estimated and found 
to dominate the rate. Therefore, in the present calculation we will take into account 
effects of the strong interactions in C'|®(mc). On the other hand, as mentioned above, the 
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operator Og mixes with four-quark operators even in the absence of QCD corrections |T^ . 
Finally, the RG running does not affect Oio, i.e. Cio{mc) = Cio{Mw). Thus, in order 
to estimate the c — ui^i~ amplitude it is a good approximation to consider the QCD 
effects only where they are dominant, i.e. in whereas we expect these to be less 

dramatic in Cf^{mc). 

The leading order mixing of with Og results in 


a 


(') eff 


= cPiMw)+ E 


i=d,s,(b) 


2 mj 

In- 


+ 


8^2 


9 9 s 


where we have dehned 






1 - 


4^2 


'Tiz,) 


( 10 ) 


T{z) = 


2 arctan 


--1 


In 




_ 4z^ 


(for s < 4^2) 


( 11 ) 


in (for s > 4^2) ^ 


and s = sjml, Zi = mi/rric. The logarithmic dependence on the internal quark mass m* 
in the second term of Eq. (|I^) cancels against a similar term in the Inami-Lim function 
Fi (xi) entering in Cg{Mw), leaving no spurious divergences in the ruj —>■ 0 limit. 


To compute the differential decay rate in terms of the Wilson coefficients, we use 
the two-loop QCD corrected value of Cj^{mc) as obtained in Ref. |P, compute C|®(mc) 
from Eq. (pR]), and Cio(mc) = Cig{Mw) from Eq. (y). The differential decay rate in the 
approximation of massless leptons is given by 


dTc. 


«£+£- 


ds 


-|-12 Cp{mc) Re 


cP ' + ICio 


G 


(') eff 


{rric) -h 4 (1 -h T 


t) (1 + 2S) 

Cf{m,) " 


( 12 ) 


where td refers to the lifetime of either or D^. We estimate the inclusive branching 
ratios for rric = 1.5 GeV, rus = 0.15 GeV, mb = 4.8 GeV and = 0, 


^ D+- 




~ 2 X 10” 


Br 


(sd) 
D0->X?.e+e 


_ ~ 8 X 10 


-9 


(13) 


It is useful to observe that the dominant contributions to the rates in Eq. ([T^) come from 
the leading order mixing of Og with the four-quark operators oP, the second term in 
Eq. (pR|). As noted above, the dominance of light-quark intermediate states in the short 
distance contributions is a signal of the presence of large long distance effects. However, 
when considering the contributions of various new physics scenarios, it should be kept in 
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Table 1: Examples oi D ^ PV^ — P^~^^ Mechanism. 


Mode 

g^(pole) 

g^(expt) 

D~^ —>■ TT^cj) —> 7r^e^e“ 

1.8 ■ 10“® 

< 5.2 ■ 10-6 


1.5 ■ 10-6 

< 1.5-10-6 

D+ —>■ 7r+0 — > 7r+e+e“ 

1.1 ■ 10-6 

< 2.7-10-^ 


0.9-10-6 

< 1.4-10-^ 


mind that their magnitudes must be compared to the mixing of these operators. Shifts in 
the matching conditions for the Wilson coefficients C 7 , Cg and Cio, even when large, are 
not enough to overwhelm the long distance effects in most extensions of the SM. These 
considerations will be helpful when we evaluate what type of new physics scenarios might 
be relevant in these decay modes. 


2.1.2 The Long Distance Contributions to D ^ 


As a hrst estimate of the contributions of long distance physics we will consider the 
resonance process D XV where V = (j),p,uj. We isolate contributions from 

this particular mechanism by integrating dT / dq^ over each resonance peak associated with 
an exchanged vector or pseudoscalar meson. The branching ratios thus obtained (we refer 
to each such branching ratio as are in the 0(10“®) range. Modes experiencing 

the largest effects are displayed in Table |l] (see also Ref. JT^), where we compare our 
theoretically derived branching ratios with existing experimental bounds . Due to the 
small r] —»• and rj' —>■ branching ratios, the dominant contributions arise from 

exchange. 

This result suggests that the long distance contributions overwhelm the short distance 
physics and possibly any new physics that might be present. However, as we will see 
below this is not always the case. A more thorough treatment requires looking at all the 
kinematically available regions in D —>■ Xui~^i~, not just the resonance region. In order 
to do this, the effect of these states can be thought of as a shift in the short distance 


coefficient in Eq. (IT), since V —> selects a vector coupling for the leptons. This 


follows from Ref. ii , which incorporates in a similar manner the resonant contributions 
to 6 —*• decays via a dispersion relation for —>■ hadrons. This procedure is 

manifestly gauge invariant. The new contribution can be written via the replacement 


/Oeff 

^9 




^£+e- 


^ — s — imyd 




Vi 


(14) 
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where the sum is over the various relevant resonances, my. and T y, are the resonance mass 
and width, and the factor Kj ~ ^iX) is a free parameter adjusted to £t the non-leptonic 
decays D XVj when the Vj are on shell. We obtain ~ 3.6, Kp ~ 0.7 and ~ 3.1. 
The last value comes from assuming = 10“^, since a direct measurement is not 

available yet. 

As a hrst example we study the > 7 r’''e’''e“ decay. The main long-distance contri¬ 
butions come from the 0, p and uj resonances. The rj and rj' effects are negligibly small. 
The dilepton mass distribution for this decay takes the form 

+ |CioP^ , (15) 


dT 

ds 


1927r5 


ip.n/+(5)p 


ZrUr 


mn 


Cf + Cf 


where s = m^g is the squared of the dilepton mass. Here we have make use of the heavy 
quark spin symmetry relations that relate the matrix elements of O 7 to the “semileptonic” 
matrix elements of Og and Oig |]^. An additional form-factor is formally still present, 
but its contribution to the decay rate is suppressed by (m^/m^))^ and is neglected here. 
For the form-factor /+(s) we make use of the prediction of Chiral Perturbation Theory 
for Heavy Hadrons [^, which at low recoil gives 


f+X) = 


fp gD*Dn 

UX-s/MX) 


(16) 


where we use the recent CLEO measurement gD*Dn = 0.59 ± 0.1 ± 0.07, and we take 
Jd = 200 MeV. In Fig. |l| we present this distribution as a function of the dilepton mass. 
The two narrow peaks are the 0 and the uj, which sit on top of the broader p. The total 
rate results in Brp+^^+^+g- ~ 2 x 10“®. Although most of this branching ratio arises 
from the intermediate 7r~^(j) state, we can see from Figure |l]that new physics effects as low 
as 10 “^ can be observed as long as such sensitivity is achieved in the regions away from 
the UJ and (j) resonances, both at low and high dilepton mass squared. 

Similarly, we can consider the decay Since there is less data available 

at the moment on the D —> VV modes, we will take the values of the Ki in Eq. (pAf) from 
the fits to the —>■ tt+H case studied above. For the semileptonic form-factors we use the 

extracted values from the D K*iiJ data and assuming SU{3) symmetry^. The total 
integrated branching ratio is Hr^io^pOg+e- = 1.8x 10“® (i.e. Brp+^p+^+e- =4.5x10“®). As 
can be seen in Fig. H, once again most of this rate comes from the resonance contributions. 
However, there is also a region -in this case conhned to low values of ruee due to the 
kinematics- where sensitive measurements could test the SM short distance structure 
of these transitions. In addition, the p modes contain angular information in the form 
of a forward-backward asymmetry for the lepton pair. Since this asymmetry arises as 


^The > /9 form-factors will be extracted with precision at Charm and B factories. In the meantime, 
we do not believe the assumption of SU (3) symmetry will affect our main conclusions here. 
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D*->-rT* e* e 



Figure 1: The dilepton mass distribution for —»• normalized to Tu+. The 

solid line shows the sum of the short and the long distance SM contributions. The dashed 
line corresponds to the short distance contribution only. The dot-dash line includes the 
allowed R-parity violating contribution from Supersymmetry (see Section |3.1.2|) 




>p° e* 


e 



Figure 2: The dilepton mass distribution for > p^e~^e~, normalized to Fjjo. The solid 
line shows the sum of the short and the long distance SM contributions. The dashed 
line corresponds to the short distance contribution only. The dot-dash line includes the 
allowed R-parity violating contribution from Supersymmetry (see Section |3.1.2|) 
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a consequence of the interference between the vector and the axial-vector couplings of 
the leptons, it is negligible in the SM since the vector couplings due to vector mesons 
overwhelm the axial-vector couplings. This is true even away from the resonance region, 
partly because of the large width of the p and partly since the coefficient cl and Cj 
get large enhancements due to mixing with O 2 and from the QCD corrections, whereas 
Cio -the axial-vector coupling- is not affected by any of these. This results in a very small 
interference. We expand on this point and consider the possibility of large asymmetries 
from physics beyond the SM in Section p.l.2| . For both the tt and p modes the sensitivity 
to new physics effects is reserved to large 0{1) enhancements since the long distance 
contributions are still important even when away from the resonances. 


We hnally compare our results in Figs. ^ and with those obtained in Refs. 
and [O. The short distance calculations in both these papers do not include the tree- 


level mixing of O 9 with O 2 . This effect determines most of the short distance amplitude. 
Also, as mentioned above, this piece cancels the logarithm in Eq. ( 0 ). a scheme dependent 
term of no physical signihcance. If this cancellation did not take place the logarithm would 
be the largest contribution to Cg. In addition, in Ref the QCD corrections are not 
included. We also differ in the long distance results, which dominate these decays. For 
D —> the authors of Ref. make use of the factorization approximation, as well 

as heavy hadron chiral perturbation theory for both pseudoscalars and vector mesons. 
It is far from clear that the use of both approximations in D decays is warranted. For 
the case of D —»• pl'^l~ ^ the results of Ref. show a large enhancement at low (f 
when compared with Fig. However, a enhancement can only appear as a result 
of non-factorizable contributions. This is clear from Ref. and |^: the factorization 
amplitude for D pV, when combined with a gauge invariant (7 — V) mixing, leads 
to a null contribution to D —>■ . This is due to the fact that the mixing of the 

operator O 2 with Oj is non-factorizable [^. A resonant contribution to 07 , leading to a 
1/q^ behavior, is then proportional to which is mostly given by the O 2 mixing. In 
addition, when compared with the usual short distance matrix element of O 7 , this resonant 
contribution will be further suppressed by the factor , where gvi^f) is the 

(7 — V) mixing form-factor, and parametrizes the non-factorizable amplitude 

{pV\Oy\D), which is of O^K^Qj^/rric) @|. Thus, even if we take the on-shell values for 
these quantities, the resonant contribution to O 7 is likely to be below 10% of the SM short 
distance contribution. The actual off-shell values at low far from the resonances are 
likely to be even smaller. We then conclude that the 1/q^ enhancement is mostly given by 
the short distance contribution. This is only noticeable at extremely small values of the 
dilepton mass, so that it is likely to be beyond the experimental sensitivity in the electron 
modes (due to Dalitz conversion), whereas in the muon modes it lies beyond the physical 
region. On the other hand, the factorizable pieces contribute to the matrix elements of 
O 9 , just as in Eq. (H), and give no enhancement at low values of q^. 
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(a) (b) 

Figure 3: Some long distance contributions. 

2.2 Neutrino-antineutrino Emission D 

In the Standard Model, decays such as 

L)+(p) ^ 7r+(p') Pi(k) and F>°(p) ^ K^{p) vip) ui{k) (17) 

will have branching ratios which are generally (but, as we shall show, not always) too 
small to measure. Such decays thus represent attractive modes for new physics searches. 


2.2.1 The Short Distance Contribution c —> uiyeue 


These decay modes are induced by Z penguin as well as box diagrams. The corresponding 
effective hamiltonian takes the form 

H {A,X^(a:,) + XbX^Xb)] ■ ( 18 ) 


The functions in Eq. (^) are defined by X^{xi) = D{xi,me)/2, with the functions D 
given in Ref. |^. Although we have explicitly kept the dependence on the charged 


lepton masses arising from the box diagrams, this is of numerical significance only when 
considering the strange quark contributions with an internal tan lepton. In any case, the 
branching ratios in the SM are unobservably small. For instance, one has 


Br 


(s.d.) 

D+^XuVu 


~ 1.2 X 10 


-15 


Br 


(s.d.) 

DO^XuUp 


~ 5.0 X 10“^® 


(19) 


where the contributions of all neutrinos have been included. 


2.2.2 Long Distance Contributions to D ^ Pr'eJ^i 

Long-distance contributions to the exclusive transition D —> PuePe {P is a pseudoscalar 
meson) can have just hadrons, just leptons or both hadrons and leptons in the intermediate 
state. Examples of the first two cases are depicted respectively in Fig. ^a) and Fig. ^(b). 

As a simple model of the purely hadronic intermediate state, we consider in detail 
the non-leptonic weak process D{p) 7r(p')I/®(g) followed by the conversion r“(«) ^ 
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ue{k)u£{k), c/Fig. |](a). We determine first the (V^ = 0,p°,a;) vertex, which 

has the invariant amplitude 


A4vo- 




92 


2 cos( 


1 


-j u{k)T^]:<-’m (oi E-z? V") . 


( 20 ) 


where is the current coupling quark q to the Z gauge boson. Only the vector part of 
the current contributes and we hnd 


Af yo_ 




2Gf 

72 


■hvu{k)eyT^v{k) 


( 21 ) 


Using the measured electromagnetic transitions U° —>■ i~^£ (U° = p^,u!,(j)) as input, we 

hnd for the coupling hy 


(3/2 - 2sl)Ml/U ^ 0.112 GeV' 


(U = 0) 


|hy| = <1 (9/8 - 2sl)Mllfp - 3M^/8U - 0.107 GeV^ (U = p) 
-(9/8 - 2sl)Mllf^ + 3MV8/, ~ 0.008 GeV" (V = uj) 


( 22 ) 


where we adopt the numerical values of f^, fp, listed in Ref. 


The corresponding transition amplitude for the non-leptonic D decay process is then 

,(V 0 ) , ,0 1 


= GlMl- 


q^-{My-iYy/2Y 


F{(:^)hy{(^)u{k)p ■-iYip}{k) , (23) 


where g=p — p' = fc + /cis the four-momentum carried by the virtual vector meson and 
F{q^) appears in the D V^P amplitude. We hnd for the g^-distribution 


G%M% Ip'I F\q^)hl{q^) 


dq^ 


1927r3 Ml (g2 - + F^M^ 


{q-p 


f\2 


q^M^ 


V 


(24) 


We have used data from non-leptonic decays into pseudoscalar-vector hnal states (D 
P + U°) to serve as input for —>• tt+z/^z/^ (p° pole), R'^z/^z/^ {p^,uj,(j) poles) and 

Dt —>• poles). Taking the largest contributor in each category, we obtain 


~ 5.1 X 10 
Br ~ 2.4 X 10 

^ 7.8x10-^^ 


(U = p°) 
(U = 0) 
{V = <P) , 


(25) 


where we have summed over the three neutrino havors. Although this analysis pertains 
to just the amplitudes of Fig. ^(a), we believe our results rehect the order of magnitude 
to be expected for other hadronic intermediate states as well. All such processes lead to 
unmeasurably small branching ratios. 
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There will also be amplitudes with single lepton intermediate states, as in Fig. ib). 
For electron and muon intermediate states, the amplitude for D{p) P{p')v(^{k)v^{k) is 
reducible to 


= -2GlV^,V:,u(k)p ■ tFlU*) + . (26) 

These lead to the branching ratios 


BrD+-.n 




~ 1.8 X 10"^® , 


Br 






~ 3.8 X 10"^® 


(27) 


which are again too small for detection. 

There remains the case in which propagates as the intermediate state. This differs 
from the above cases involving e and p propagation in that for part of the Ut-Vt phase 
space, the intermediate is on the mass shell. The mode —> r’*' + has been 

observed^ with = (7 ± 4)% whereas ^ + i/,- has not (the predicted 

branching ratio is ~ 9.2 10“^). Once the on-shell has been produced, its 

branching ratio to decay into a given meson can be appreciable, e.g. Br^-^p+p^ ~ 0.25, 
Br-r^TT+u-r — 0.11, etc. Such transitions, although involving production of a uu pair in the 
hnal state, should be measurable at a 5 and/or Charm factory. 


2.3 Two Photon Emission 77 

The amplitude for the transition D^{p) —> 7 (^ 1 , Ai) 7 (g 2 , A 2 ) can be expressed as 


= et (l)et (2) 



qi ■ q2 g^'') Cdo^^ + 


(28) 


The invariant amplitudes i?_D 077 and Cno^^ are P-conserving and P-violating, respectively, 
and carry units of inverse energy. They contribute to the —> 77 branching ratio as 


Br £)0^^.y — 


- 


I-8^077 P + ICDO77/ 


Odvr 

The amplitude in Eq. (^) is sometimes written in the equivalent form 


M 


D^'y'y 


Co^ry-y 


FrF2^u+i 


B 


Z)®77 


2 1 I 2 

where F^'^ = — q'^e^) and 

^In this experiment, only the leptonic decay mode r'*" - 


F\puF2p,u 1 


sVt was detected. [Ea 


(29) 


(30) 
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(a) (b) 

Figure 4: IPi? contributions to c —^ M 77 . 


2.3.1 The Short Distance Contribution cu 


77 


Consider the quark level transition c —M 77 . This can arise via one-particle irreducible 
(IPI) processes in which both photons arise from the interaction vertex or one-particle 
reducible (IPR) processes in which at least one of the photons is radiated from the initial 
state c-quark or hnal state n-quark. 

To estimate the c —> M 77 amplitude, we employ an approximation which makes use of 
known results on the related process c —>• nq. According to Ref. |^, the two-loop c —>■ wy 
vertex is 

, (31) 


w(s.d.) ^ 

y/2 167r2' 


where |A| ~ 0.0047. Keeping in mind that there are additional diagrams which must 
be accounted for in a complete two-loop analysis, we shall use this as input to the IPR 
graphs depicted in Fig. ^ The dominant contribution to the c —> M 77 amplitude involves 
photon emission from the n-quark. To ensure that the effect is indeed ’short-range’, we 
follow the locality procedure employed in Ref. . This yields for cu 77 the amplitude 


|^(s.d.) I _ |^(s.d.) I _ 






G pcx 


mr_ 


D^'yy I 


Md - rric 


fD |A| 


resulting in the branching ratio 


- 3 X 10"“ 


(32) 


(33) 


for the choice Mp — rric — 0.3 GeV. 


2.3.2 Long Distance Contributions to 


77 


We shall model long-distance contributions to the —>■ 77 amplitude using the vector 
meson dominance (VMD) mechanism and the unitarity constraint. The latter can only 
be done in a limited context since there will be many unitarity contributions. We will 
consider several one-particle intermediate states (as used in iP —77 decays) as well as 
the two-particle K^K~ intermediate state. 
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Figure 5: Vector dominance (VMD) contribution. 


Vector Meson Dominance 


One can view {c.f. Fig. |^) the —>■ 77 amplitude as the single VMD process 

+ E Vk* ^ 7 + 7 . (34) 

k 

We have previously used the VMD mechanism to model the general single-photon emission 
D —>■ M-I-7 (M is some non-charm meson) |^. It is straightforward to extend our analysis 
to the —> 77 mode, as long as care is taken in the 77 amplitude to ensure gauge 

invariance and Bose-Einstein statistics. The amplitudes used in the D^{p) —>■ V^{k)+'j{q) 
transition are defined as 


■MdV't ~ ^ 7 ) 


Cv{.Kq^l 


k ■ qSia,) + 


(35) 


The VMD amplitude that we calculate is therefore of the form 


p(vmd) 

^DO'yy 


^ 2e 

E 7- Bvi Vi , 

i JVi 


c 


(vmd) 




r 

-r~ 

JVi 


Vi 


(36) 


where fv is the coupling for the — ■y conversion amplitude, the index 7’ refers to the 
specific vector meson 0°) and 7 * is a factor accounting for the VMD extrapolation 

made in We take 7 ^ ~ 1/2 as a reasonable choice. 

The values in Table are somewhat lower than those which would be obtained from 
the V 7 amplitudes in Ref. [^]. The main reason for this is the central value for BrDo^^pO, 
which is a numerically signihcant input to the VMD calculation, cited in the Particle 
Data Group compilation has decreased by a factor of about three between 1994 and 2000. 
Using the central values in Table |] and assuming positive interference between the various 
amplitudes to provide the maximal VMD signal gives the branching ratio 

= (3.5 i 7 ) X10-* . ( 37 ) 
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Table 2: VMD Amplitudes (10“® GeV-^). 



^ l/°7 

Dvmd 

■^D^77 

/^vmd 

po 

^ P°7 

0.036 (1 ±0.7) 

0.045 (1 ±0.3) 


—>■ a;°7 

0.011 (1 ±0.5) 

0.012 (1 ±0.5) 


^ 0% 

0.047 (1 ±0.7) 

0.036 (1 ±0.4) 


Single-particle Unitarity Contribution 


In this category of amplitudes (c/. Fig. P) the mixes with a spinless meson (either 
a pseudoscalar or a scalar 5„) and finally decays into a photon pair, 


B 


(mix) 


1 


DO,, = 

Pn "'"P* 

(mix) _ /c I'iv(p-V-)| nO' ^ 


Ml - Mp2 


/c i'X7''P''''h n0\ _t 

^D0„ Z. W ) Ml-Ms2 ■ 

On ^ ^ri 

Let us consider two distinct kinds of contributions, 

1. If the spinless meson is a ground-state particle (tt^, r] or ?7'))EI have 

)gnd ^F^2fD^ 


(38) 


B 


D®77 




TT 


Mi 


Mj ^ 

v/2M|,-M2 372 


fk{e) 


(39) 


where 02 ^ —0.55, 9 ~ —20°, fr){9) = cos^ 0 — 2-\/2sin6'cos6' and fn'{9) = sin^^ -|- 
2v^sin0 cos6*. The above parameterization for the two-photon vertices agrees with 
the values determined experimentally. 


0.0249 GeV^ (7r°) 
Bp^,, = { 0.0275 GeV^ (r/) 
0.0334 GeV^ (7) 


(40) 


B^o,, is seen to vanish, as it must, in the limit of SU(3) flavor symmetry (there 
= 0 and the 7r°, p contributions cancel). From Eq. (^9]), we obtain 
the branching ratio 

Br^nL,, - 3 x lO-^' . (41) 


^The kaon intermediate state is disfavored due to the small K 77 branching ratio. 
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Figure 6: Weak mixing contribution. 


2. If the intermediate meson is a spinless resonance the decay chain becomes 

—>■ yy. Since little is yet known about meson excitations, both the weak mixing 
amplitudes and the two-photon emission amplitudes must be modeled theoretically. 
The D^-to-resonance weak matrix element will depend upon the flavor structure of 

e.g. 

= {uu — dd) / y/2) 

= ss) (42) 




GF(^2fD 


^dfu/ 

^sJr 

v:,VusfH 


where the flavor content of is in parentheses and estimates for resonance decay 
constants Jr are given in Ref. |^. The R^ —> yy mode has been observed for 
a number of resonances and has typical branching ratios = 0(10“®) for 

Mr ~ 1 —> 1.3 GeV, decreasing to Btrq^^^ = 0(10“®) for Mr > 1.5 GeV. 

For a concrete example of the resonance mechanism, we choose R® = 7r(1800) and 
assume i3r,r(i800)^77 — 10“®. The resulting O® —> yy branching ratio is 


„ R°=R1S00) 


10 “io 


(43) 


Two-particle Unitarity Contribution 


In a factorization approach, the O® —^ K^K amplitude {cf. Fig. is 


MnoK-^K- = 


V V* f 

''cs ''usJ 


Mi. 


U(Ml) + 


(44) 


where /± are form factors and / is a constant containing information about QGD cor¬ 
rections and the kaon decay constant. A £t to the measured —> K^K~ decay rate 

yields 


/ 



f+i^x) + 


Ml 

Ml 




141 MeV 


(45) 
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(a) (b) 

Figure 7: Unitarity contributions: (a) K~^K~, (b) K*^K*~. 

Similar to the Bg system [^] , the K^K~ intermediate state contributes via unitarity to 
only the amplitude of Eq. (^) and is proportional to precisely the same combination 

of form factors appearing in Eq. (^), 

Ml 


Im C 


from which we obtain 


{K+K-) 

D^-yy 


= 2a 


'■K 

Mh 


Br 


{K+K-) 


Summary of 77° 


1 — MIdok+k- 


0.7 X 10" 


(46) 

(47) 


77 


Considered together, the above examples lead us to anticipate a branching ratio in the 
neighborhood of 10“®. Our maximal (he. constructive interference) VMD signal has a 
central value ~ 3.5 x 10“®. The recent work of Ref. ||^ provides an independent 

estimate of the 77° —> 77 transition and obtains a similar order-of-magnitude result. 


2.4 Lepton-antilepton Emission ^ £+£- 

The general form for the amplitude describing 77°(p) —> s+)7'“(A;_, s_) is 

= u{k-, sJ)[Aj:,oi+i- + 75 77 £) 0 £+£-] n(fc+, s+) , 
and the associated decay rate is 

Md 


^ - 


Svr \ 


1-4 


m; 


I +1—4 


rrip 


I B]:,0i+i -1 


(48) 


(49) 


2.4.1 Short Distance Contributions cu —> 


The short distance {0{as) corrected) transition amplitude is given by [p!4 

GlM^fome 


ri(s.d.) 


TT^ 




(50) 
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(a) 


(b) 


Figure 8 : Unitarity contributions: (a) One-particle, (b) Two-particle 77 . 


where 


F= E 

i=d,s,b 


V v* 

''ui Cl 


Xi , as 


b 2 4: + 71^ 

In Xi H--— 


(51) 


with Xi = The amplitude A£) 0 ^+^- vanishes due to the equations of motion. 

The explicit dependence on lepton mass in the decay amplitude overwhelmingly favors 
the hnal state over that of Upon employing the quark mass values ~ 

0.01 GeV, TUs — 0.12 GeV, mt, ~ 5.1 GeV, the Wolfenstein GKM parameters A ~ 0.22, 
A ~ 0.82, p ~ 0 . 21 , T] ~ 0.35 and the decay constant /d — 0.2 GeV, we obtain the 
branching fraction — 10 “^®. 


2.4.2 Long Distance Contribntions to DO ^ 

In the following, we consider two long distance unitarity contributions (c/. Fig. H) which 
lead to —>• transitions. In each case, the decay amplitude is dependent on 

the lepton mass, and thus we shall provide numerical branching ratios only for the case 
DO ^ p+p-. 


Single-particle Unitarity Contribntion 


The single-particle ‘weak-mixing’ contribution to can be estimated in a 

manner like that considered for the DO —> 77 transition (c/. Eq. ®). For dehniteness, 
we consider the D'^ —^ parity-conserving amplitude D£) 0 £+£- (see Eq. 


n(mix) 

^D0l+£- 


Pn 


1 

Ml - Mp. 




(52) 


and we write ground state {n^,r],p') and resonance 

contributions. 

There is little known regarding the PnP~^p~ {Pn = t^o, 7 , p') vertices. In the following, 
we assume these quantities have the same flavor structure as the corresponding Dn 77 ver¬ 
tices described earlierand obtain the overall PnP^p~ normalization from the measured 
^This ensures that our expression will vanish in the limit of SU(3) flavor symmetry. 
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T] /i"*"/!" mode. From this we predict for the r]' (960) —>■ mode a branching ratio 

~ 5.6 X 10“^, well below the cnrrent bonnd Br^i^+^- < 10“^. The gronnd state 
contribntion is then 


B 


(gnd) 

D0£+£- 


GFO‘2fDBp 




72 

26 - 6 


^yj-rj 


372 M2 -M2 
26 - 6 


6 M 2 
L 72 M2 - M2 


(^cos^ 6 — 2A/2sin6* cos^ 


M'2 / , r- 

^ , (sin^ 6 + 272 sin 6 cos 6 

— A/f ^ V 


372 Ml-M‘ 

with Bp^+^- = 3.47 x 10“®. This leads to the branching ratio 


Br 


(gnd) 

D°^£+£- 


~ 2.5 X 10 


-18 


(53) 


(54) 


There can also, in principle, be intermediate state contribntions from = 0^ nentral 


resonances {7?°}. Using the D'^-io-R^ mixing amplitnde already obtained in Eq. 
again identifying the resonance as 7r(1800), we hnd 


- 1-8 X 10 


3-*- 7r(1800)f+£- 


M. 


= 1.8 X 10“^i3r 


(1800) 


'7r(1800)^r+£- 


Upon assnming i3r^(i8oo)^£+£- = 10 ^2 as onr defanlt branching ratio, we obtain 


».(^(1800)) _ c: n V in-17 '^^^(1800)^£+r 

— o.u X iu 


and 


(55) 


( 66 ) 


Althongh possibly enhanced relative to the light-meson pole contribntions, the resnlt is 
still nnmeasnreably small. 

The Two-photon Unitarity Contribution 


In the Kl ^ e^e transition, the two-photon intermediate state is known to play an 
important role. Let ns therefore consider the contribntion of this intermediate state for 

£+r. 


Xm 

Ai,A 2 ' 


(Rqx (Rq2 


2a;i(27r)3 2uj2{2T£f 


X M*^^f,+^-{27r)^6^‘^\p -qi- ^ 2 ) 


Upon inserting the general form of the 77 appearing in Eq. (|^ , we obtain 


Im A^2o}+f- = arriiB 


Ml 


m 


2 ’ 


A(f 2 

Xm = iamiCpo^^ In —y 

Tfln 


We hnd 


^ 2.7 X 10-^i3r,,o_ 


>77 


(67) 

(58) 

(59) 

(60) 
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Summary of 

The largest of our estimates, the two-photon unitarity component, for the long distance 
contribution to jj,'^fx~ favors a branching ratio somewhere in excess of 10“^^. More 

generally, it scales as 2.7 x 10“® times the branching ratio for 77 . With the estimate 

> 10 “® arrived at in the previous section, we therefore anticipate a branching 
ratio for of at least 3 x 10“^^. 

3 Potential for New Physics Contributions 

As discussed in the introduction, the charm system provides a unique laboratory to probe 
physics beyond the Standard Model as it offers a complementary probe of physics to that 
attainable from a study of rare processes in the down-quark sector. As we found in the 
previous section, short distance SM contributions to rare charm decays are quite small 
due to the effectiveness of the GIM mechanism, and most reactions are dominated by 
long range effects. However, we saw that for some reactions there exists a window for 
the potential observation of new short distance effects, in particular for specific regions 
of the invariant dilepton mass spectrum in D ^ Indeed in some cases, it is 

precisely because the SM rates are so small that charm provides an untapped opportunity 
to discover new effects and offers a detailed test of the SM in the up-quark sector. 

In this section, we delineate some new physics possibilities, motivated by supersym¬ 
metric, grand-unihed, extra dimensional, or strongly coupled extensions of the SM, which 
give rise to observable effects in rare charm transitions. In some cases, we hnd that present 
experimental limits on these channels already constrain the model parameter space. 

3.1 Supersymmetry and Rare Charm Decays 

We hrst examine the effects of Snpersymmetry (SUSY) in rare charm decays, concentrat¬ 
ing on the exclusive modes D and D —> Weak scale Snpersymmetry 

is a possible solntion to the hierarchy problem and as snch is a well motivated theory 
of physics beyond the SM. We consider the general case of the unconstrained version 
of the Minimal Supersymmetric extension of the Standard Model where no particular 
SUSY breaking mechanism is assumed and investigate the two scenarios where R-parity 
is conserved or violated. Imposing the constraints on the SUSY parameter space from 
cnrrent data, we hnd that in both cases, the snpersymmetric contribntions to these decay 
channels can be qnite large, particularly in the low dilepton mass region [i.e. below nip). 


20 



c 



Figure 9: A typical contribution to c ^ u FCNC transitions in the MSSM. The cross 
denotes one mass insertion (^i 2 )aa') with A, A' = L, R. 


3.1.1 Minimal Supersymmetric Standard Model 


The Minimal Supersymmetric Standard Model (MSSM) is the simplest supersymmetric 
extension of the SM and involves a doubling of the particle spectrum by putting all SM 
fermions in chiral supermultiplets, as well as the SM gauge bosons in vector supermul- 
tiplets. In our discussion, we do not assume any particular Supersymmetry breaking 
mechanism, but rather use a parameterization of all possible soft SUSY breaking terms. 
A large number, of order 100, of new parameters is then introduced. The soft super- 
symmetry breaking sector generally includes three gaugino masses, as well as trilinear 
scalar interactions, Higgs and sfermion masses. Supersymmetry contains many potential 
sources for flavor violation. In particular, if we choose to rotate the squark helds by the 
same matrices that diagonalize the quark mass matrices, then the squark mass matrices 
are not diagonal. In this super-CKM basis, squark propagators can be expanded so that 
non-diagonal mass terms result in mass insertions that change the squark flavor |]^,|30 
These mass insertions can be parameterized in a model independent fashion via 




M? ’ 

±y±q 


(61) 


where i ^ j are generation indices. A, A' denote the chirality, {M^jY are the off-diagonal 
elements of the up-type squark mass matrix, and Mq represents the average squark mass. 
The exchange of squarks in loops thus leads to FCNC through diagrams such as the one 
depicted in Fig. |^. This source of flavor violation can be avoided in specihc SUSY breaking 
scenarios such as gauge-mediation or anomaly mediation, but is present in general. It 
appears, for instance if SUSY breaking is mediated by gravity. 

The MSSM contributions to loop mediated processes in addition to those of the SM 
are: gluino-squark exchange, chargino/neutralino-squark exchange and charged Higgs- 
quark exchange. This last contribution carries the same CKM structure as in the SM 
loop diagram and is proportional to the internal and external quark masses; it thus leads 
to small effects in rare charm transitions and we neglect it here. The gluino-squark contri¬ 
bution proceeds via flavor diagonal vertices proportional to the strong coupling constant 
and in principle dominates the CKM suppressed, weak-scale strength chargino/neutralino- 
squark contributions. We thus only consider the case of gluino-squark exchange here as 
an estimate of the potential size of supersymmetric effects in rare charm decays. We note 
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that the analogous gluino contributions to rare K and B transitions have led to strong 
universality constraints on the charged Q = —1/3 squark sector [^]. Here, we examine 
the level at which the corresponding constraints can be obtained in the charged Q = +2/3 
squark sector once data accumulates at B and charm factories. 

Within the context of the mass insertion approximation the effects are included in the 
Wilson coefficients corresponding to the decay D —>■ X£^£~ via 


Q = cr + Cf , 


(62) 


for z = 7, 9,10. Allowing for only one insertion, the explicit contributions from the gluino- 
squark diagrams are [E 


= “9 


£ XU \ -^132 (^) \ £ XU \ u I 


rrir 


(63) 


and 


Cl = - 


V2 


27 GpMl 




(64) 


with the contribution to Cio vanishing at this order due to the helicity structure. If we 
allow for two mass insertions, there is a contribution to Cg lo given by 


1 O' 

Cfo ^ ~Q~ {S22)LR{5u)LRPo32iu) = 

9 a 


c. 


1-4 sim d 


w 


(65) 


Here, u = M|/M| and the functions Pijkiu) are dehned as 


Pijk{u) = dx 


x\l — xY 


0 (1 — a: + uxY 


( 66 ) 


In addition, the operator basis can be extended by the “wrong chirality” operators O 7 , O 9 
and Oio, obtained by switching the quark chiralities in Eqs. (|) and (|^. The gluino-squark 
contributions to the corresponding Wilson coefficients are 


C'l = -- 


Cl = 


V 2 

9 CpM? 

V2 


TTOo 


(5: 


Pl32iu) , 
12 )rR -^ 


U 


12)LRJP122 


P22(U)^ 
rrir 


27 GpMl 


nas ((5i2)rr-Po42(m) - (1 - 4sin^ ^w)C'io , 


Ag _ 

'-^10 — “ 


1 ^ 
9 a 


{^22)lr{^12)lrPo32{u) 


(67) 


where the expression for (Tfo is again obtained with a double insertion. 

As was noted in Refs. |3^J33[| , in both Cf and Cf the term in which the squark chirality 
labels are mixed introduces the enhancement factor Mg/iric- In the SM the chirality flip 
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which appears in O 7 occurs by a flip of one external quark line, resulting in a factor 
of rric included in the operator’s dehnitionQ. However, in the gluino-squark diagram, the 
insertion of (^i 2 )hl forces the chirality flip to take place in the gluino line, thus introducing 
a Mg factor instead of rric. This yields a signihcant enhancement in the short distance 
contributions to the process D Xu'y |] 3 ^, which is unfortunately obscured by the large 
long range effects. 

The most stringent bounds that apply to the non-universal soft breaking terms (^i 2 )aa' 
come from the experimental searches for mixing. The current CLEO limit 

implies 


/Am 


D 


cos 5 -|- 


f—y 

\ 2 r £)0 ) 


sin(5 > < 0.04% 


( 68 ) 


2 [ V T/jo 

where 5 is a strong relative phase between the Cabibbo-allowed and the doubly Cabibbo- 
suppressed Kn decays. Neglecting this phase results in the constraints obtained 

in Ref. |^, which we collect in Table These bounds were obtained assuming that 


Mi/M! 


(*^ 12 )^^ 

0.3 

0.03 

0.04 

1.0 

0.06 

0.02 

4.0 

0.14 

0.02 


Table 3: Bounds on (^^ 2 )^^; (^ 1 ^ 2 )^^ from mixing |Q (neglecting the strong 

phase). All constraints should be multiplied by (M^/SOO GeV). 


{.^12)RR — 0 (^1^2)^^? = ('^i2)hl; these assumptions are found to be numerically unim¬ 

portant. 

In order to estimate the effects in c — ui~^£~ transitions from the gluino contributions, 
we need to specify Mg and Mg. We consider four sample cases: (I): Mg = Mg = 250 GeV; 
(II): Mg = 2Mg = 500 GeV; (III): Mg = Mg = 1000 GeV and (IV): Mg = (1/2) = 

250 GeV. We hrst examine D~^ —> 7r''“e“'“e“. In Fig. 0 we show the dilepton mass distri¬ 
bution as a function of the dilepton mass. Although the net effect is relatively small in 
the integrated rate (an increase ~ 20% or smaller), the enhancement due to the SUSY 
contributions is most conspicuous away from the vector resonances, particularly for low 
dilepton masses. Experiments sensitive to the dilepton mass distribution at the level of 
10“^ — 10“® can detect these SUSY contributions. However, the decays to a vector me¬ 
son, such as D pe'^e~, are more sensitive to the gluino exchange, as can be seen from 

^ The niu term, proportional to the (1 — 75 ) in the operator, is neglected. 

^Limits obtained from charge and color breaking (CCB) and bounding the potential from below 
(UFB) apply to the trilinear terms but not to the squark mass terms. Thus, unless the squark mass 
matrices are kept diagonal, CCB and UFB arguments cannot be used to constrain the non-universal mass 
insertions. 
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Figure 10: The dilepton mass distribution for ^ 7 r’''e’''e (normalized to r£)+), in 
the MSSM with non-universal soft breaking effects. The solid line is the SM; (I): Mg = 
Mg = 250 GeV; (II): Mg = 2 Mg = 500 GeV; (III): Mg = Mg = 1000 GeV and (IV): 
Mg = ( 1 / 2 ) Mg = 250 GeV. 


Fig. Q. The effect is quite pronounced and almost entirely lies in the low ruee region. 
This is mostly due to the contributions of (^ 12 )^^ to ^7 Eqs. ( |63D and (|67D, 

which contain the Mg/rric enhancement as discussed above. This effect is intensified at 
low = rrigg due to the photon propagator (see for instance Eq. (P^) for the inclusive 
decays). This low enhancement of the O 7 contribution is present in exclusive modes 
with vector mesons such as D but not in modes with pseudoscalars, such as 

D —>■ since gauge invariance forces a cancellation of the 1/q^ factor in the latter 

case {e.g., see Eq. (^)). This is apparent from a comparison of the low dilepton mass 
regions in Figs. (0) and (ITT]). 


We conclude that the D —> decays are considerably sensitive to non-universal 

soft breaking in the MSSM. The largest effect is obtained in case (IV) (dashed line in 
Fig. (PTD) and yields BvDo^pOe+e- — 1-3 x 10“®, which is roughly a factor of five times 
larger than the SM prediction given in Sect. | 2 . 1 . 2 | . The current experimental bound on 
this channel is 'BT’^^o/^pOe+e- < ^ 10“^. For muon final states, the somewhat more 


stringent constraint < 2.2 x 10 ® should be compared to Br^o. 


1.3 X 10 ® obtained in case (IV). Thus, searches for rare charm decays with sensitivities 
of 10“® and better will soon constrain the MSSM parameter space or observe an effect. 
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Figure 11: The dilepton mass distribution for —>■ p^e^e (normalized to F^jo), in 

the MSSM with non-universal soft breaking effects. The solid line is the SM; (I): Mg = 
Mg = 250 GeV; (II): Mg = 2 Mg = 500 GeV; (III): M~g = Mg = 1000 GeV and (IV): 
Mg = ( 1 / 2 ) Mg = 250 GeV. 


3.1.2 R Parity Violation 


The assumption of i?-parity conservation in the MSSM prohibits baryon and lepton num¬ 
ber violating terms in the super-potential. However, other symmetries can be invoked to 
prohibit rapid proton decay, such as baryon-parity or lepton-parity [^], and hence allow 
for R parity violation. The R-parity violating super-potential can be written a^ 


Wr, = \ + KikLtQ’iDk + Dl 


(69) 


where L, Q, E, U and D are the chiral super-helds in the MSSM. The SU{3) color indices 
are denoted by a ,/?,7 = 1,2,3, the SU{2)l indices by a, 6 = 1,2 and the generation 
indices are i,j, fc = 1,2,3. The helds in Eq. ( |69D are in the weak basis. The term 
is the one which is relevant for the rare charm decays we consider here as it can give 
rise to tree-level contributions through the exchange of squarks to decay channels such 
as D ^ X£~^£~, D —> as well as the lepton-flavor violating D —>■ Xp~^e~ and 

D —> p~^e~ modes. Before considering the FGNG effects in D decays, we need to rotate 
the helds to the mass basis. This leads to 


Wh, = X^jk [NiVjiDi - E,Uj] Dk + --- 

^We ignore bilinear terms which are not relevant to our discussion of FCNC effects. 


(70) 
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where V is the CKM matrix and we dehne 


. ( 71 ) 

Here, and are the matrices nsed to rotate the left-handed np- and right-handed 
down-qnark helds to the mass basis. Written in terms of component helds, this interaction 
now reads 


wv = A7qv',,|i>i.44 + <?L44 + (4)‘(4)'4l 
“444 — 444 — ( 4 ) ( 4 )'^ 4 } ■ 


(72) 


The last term in Eq. 0 can give rise to the processes c —>■ u££^'^ at tree level via the 
exchange of a down-sqnark. This leads to effects that are proportional to \'i 2 kKik with 
i = 1,2 (dne to kinematical restrictions). 

Constraints on these coefficients have been derived in the literatnre [0. For instance, 
tight bonnds are obtained in Ref. from vr+z/P by assnming that only one R- 

parity violating conpling satisfies ^ 0. We npdate this bonnd by nsing the latest 


experimental resnlt Etk+^-k+uu = (l-STlgisi) ^ 10 which yields < 0.005. 

where only one 


However, this bonnd can be avoided in the single conpling scheme 
R-parity violating conpling is taken to be non-zero in the weak basis. In this case, it is 
possible that flavor rotations may restrict the R-parity breaking indnced flavor violation 
to be present in either the charge —1/3 or -1-2/3 qnark sectors, bnt not both. Then large 
effects are possible in the np sector for observables snch as D^-D^ mixing and rare decays 
withont affecting the down-qnark sector. In Ref. a rather loose constraint on the R- 
parity breaking conplings is obtained from mixing, which conld resnlt in large effects 
in c — u££^''' decays. Here, we will take a conservative approach and make nse of more 
model-independent bonnds. The constraints on the R-parity breaking conplings for the 
processes of interest here are collected in Table ^ from Ref. . The charged cnrrent 
nniversality bonnds assnme three generations. The vr decay constraint is given by the 
qnantity The limits obtained from D —>■ K£u were first obtained in 

Ref. 0. 



\/ 

'^I2k 


\/ 

^22k 

0.02('^) 

0.04('^) 

0.06('’) 

0.2l(") 


Table 4: Most stringent (2 ct) bonnds for the R-parity violation conplings entering in rare 
D decays, from (a) charged cnrrent nniversality, (b) Rj, and (c) D —> K£h'. See Ref. 
for details. All nnmbers shonld be mnltiplied by {rn^k /lOO GeV). 
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We first consider the contributions to c 
squarks results in the effective interaction 


. The tree level exchange of down 




K2kKlk 


which after Fierzing gives 






K2tKik . 


(73) 


4 


(74) 


This corresponds to contributions to the Wilson coefficients Cg and Ciq at the high energy- 
scale given by 

Mw 


sca = -6a. = 


sin^ e 


w 


2a2 




A'o.A' 




(75) 


If we now specify £ = e and use the bounds from Table ^ we arrive at the constraint 


6Q = < 1.10 


'A( 


12k 


'A( 


life 


0.04/ \0.02 


(76) 


Notice that these are independent of the squark mass, which cancels. Taking this upper 
limit on the Wilson coefficients results in the dot-dashed lines of Figs. |I] and ^ correspond¬ 
ing to 7r~^e~^e~ and —>■ p^e~^e~, respectively. The effect in these rates is small, 

of order 10% at most, whereas the experimental bounds are a factor of 20 above this level 
in the best case (given by the pion mode). 

On the other hand, for i = p we obtain 


SC^ = < 17.4 


'A( 


22k 


'A( 


21k 


0.21 


0.06 


(77) 


These upper limits already saturate the experimental bounds of <1.5x10 ^ 

and ^ao^oOu+u- < 2.2 x 10“^ from Refs. i35|j4l[|! Thus we derive the following new 


constraint on the product of R-parity violating couplings, 

^22k ^21k < 0.004 , 


(78) 


which arises from the D~^ —>■ p mode. This allows for potentially large effects in 


both the p and tt channels as is illustrated in Figs. and |T3 


In Figure we display the dimuon mass distribution as a function of the dimuon mass 
for > ap'^p~. The solid line, corresponding to the SM prediction and including 

both the short and long distance pieces, is clearly dominated by the latter through the 
presence of the vector meson resonances as discussed above. The dashed line includes the 
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D*->rr* /j* iU 



Figure 12: The dilepton mass distribution for D'^ normalized to Tij+. The 

solid line shows the sum of the short and the long distance SM contributions. The dashed 
line includes the allowed R-parity violating contribution from Supersymmetry (see text 
for details). 


□°->p° p* p' 



Figure 13: The dilepton mass distribution for normalized to F^io. The 

solid line shows the sum of the short and the long distance SM contributions. The dashed 
line includes the allowed R-parity violating contribution from Supersymmetry (see text 
for details). 

28 








contribution of R parity violation, taking the R-parity violating coefficients to saturate 
the above bound in Eq. ([78|). It can be seen that away from the resonances there is an 
important window for the discovery of R parity violation in SUSY theories. The situation 
is similar in the distribution, shown in Figure |^. Here, the dashed line is 

again obtained by making use of the bound in Eq. (^8|). This results in an upper bound 

8.7 X 10“®, which is still below 

the experimental limit |^| Br'^^Q_ 




Rf 

)0_>p0^+^- 

2.2 X 10-^ 


In addition to the dilepton mass distribution, this decay mode also contains angular 
information as discussed in the previous section. For instance, we can dehne the forward- 
backward asymmetry for leptons as 


fi Jfl-dx — f° -^^dx 

A / 2 \ "^0 dxdq^ -1 — 1 dxdq^ 

^FB[q ) - -- 

dxp 


(79) 


where x = cos 6*, with 9 being the angle between the and the decaying D meson in 
the rest frame. Expressions for the angular distribution dT / dxdq^ can be found in 
Ref. m for the inclusive case and in Ref. for the exclusive modes. In the SM, 


Apsiq^) ffi is negligibly small throughout the kinematic region. The reason 


for this can be seen by inspecting the numerator of the asymmetry 


^Fsiq^) ~ 4 mo k Cio 


Cf 


rUc 

q2 


g f + 


ifG 



(80) 


where k is the vector meson three-momentum in the D rest frame, and /, g, F and G are 
various form-factors. Since the SM amplitude is dominated by the long distance vector 
intermediate states, we have S> Ciq. New physics contributions that make Cio — 
will hence generate a sizable asymmetry. This is illustrated in the case at hand of R parity 
violating supersymmetry. For instance, again setting the coupling to the values given in 
Eq. ([78|) , we present the forward-backward asymmetry for —> p^p^p~ in Figure |T^. In 

order to compute the asymmetry, we make use of —>• K*iu form-factors, together with 

SU{3) symmetry and heavy quark spin symmetry^. This gives a bound on the integrated 
asymmetry of ~ 0.15. For —>■ p^e'^e~, we get ~ 0.08. Supersymmetry could 

thus produce very sizable asymmetries. In general, any non-zero value of ApB{q^) that is 
measured should be interpreted as arising from new physics. 




The effective interactions of Eq. (^) also lead to a contribution to the two body decay 
—>■ p^p~. The R parity violating contribution to the branching ratio then reads 


» ftp 


■x+p- = fh ml mo 


\ 


4m2 (A 2 


\/ 

22k^21k 


m 


D 


Odvrmj 

dk 


®See the first reference cited in Ref. [liSl. 


(81) 


29 






















Figure 14: The lepton forward-backward asymmetry for , for the bound of 

Eq. (|^). (see text for details) 


Applying the bound in Eq. (|78|) gives the constraint 





< 3.5 X 10"® 




(82) 


The current experimental limit BrD 0 ^^+^- < 5.2 x 10“® is just above this value, 
implying that future measurements of this decay mode will constrain the product of these 
R parity violating couplings. 

Finally, we consider the products of R parity violating couplings that lead to lepton 
flavor violation. For instance, the products \iik^ 22 k X 2 ik^i 2 k rise to 

. This leads to 


.cr ^ ^ 4.6 X { (ll) (ll) ^ (ll) (ll) ) , (83) 

ra 

which results in < 3 x 10“®, to be contrasted with 

3.4 X 10“®. Here again, experiment is on the verge of being sensitive to R parity violating 
effects in supersymmetry. Similarly, for the corresponding two body decay we have 






< 


< 0-5 X 10-® X 


o.o2;vo.2i; vo.oe; vo.o 4 y j ’ 


(84) 


whereas the current bound is 13511 Br‘^^o_ < 8.1 x 10 ®. 


this section in Table ^ 


We summarize the results of 
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Finally, we point out that similar effects to those considered in this section are gen¬ 
erated by leptoquarks. Their exchange lead in general to effective interactions similar to 
the A' terms in Eq. (|69|). 


Decay Mode 

SM 

jRp 

Expt. Limit 


2.0 X lO-o 

2.3 X lO-o 

5.2 X lO-o 

£)0 p^e^e~ 

1.8 X lO-o 

5.1 X lO-o 

1.0 X 10-1 

— >• fX~ 

1.9 X lO-o 

1.5 X lO-o 

1.5 X lO-o 

pV+F" 

1.8 X lO-o 

8.7 X lO-o 

2.3 X 10-1 

DO ^ /i+p- 

3.0 X 10-10 

3.5 X lO-o 

4.1 X lO-o 

DO ^ e+e- 

O 

1 

to 

CO 

1.0 X 10-10 

6.2 X lO-o 

DO ^ /i+e- 

0 

1.0 X lO-o 

8.1 X lO-o 

D+ — >• 7r+p+e“ 

0 

3.0 X lO-o 

3.4 X lO-o 

DO — >■ pO^+g- 

0 

1.4 X lO-o 

4.9 X lO-o 


Table 5: Comparison of various decay modes between the SM and R parity violation. 
The third column shows how large the R parity violating effect can be. The experimental 


limits are from Refs. 16 


3.2 Extensions of Standard Model with Extra Higgses, Gauge 
Bosons, Fermions, or Dimensions 

In this section we summarize the results from classes of models which have additional 
Higgs scalar doublets, or family gauge symmetry or extra leptons. All of these give rise 
to flavor changing couplings at tree level and hence yield potentially large rates for rare 
decay modes of D mesons. In addition we briefly discuss the effects of extra dimensional 
physics on rare charm transitions. 


3.2.1 Multiple Higgs Doublets 


Many extensions of the Standard Model contain more than one Higgs scalar doublet. As 
is well known, this leads in general to tree level FCNC couplings and thus decays such 
as e+e“,/i^e^, etc may proceed at rates larger than SM expectations. In 

the down quark sector, there are severe constraints on such couplings from kaon decay 
modes |^. This does not necessarily lead to equally strong constraints on the up-quark 
sector. For example, as was shown long ago |^, it is possible that simple symmetries 
forbid AS = 1 FCNC without affecting the AC = 1 sector. 
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Let us write the general effective AC* = 1 interaction as 


Gf ^ 

ti 75 C fi(a + 675)^2 , 


(85) 


where is a model dependent dimensionless number, a and b refer to generic scalar and 
pseudoscalar couplings, respectively, and £ 1,^2 refer to the pairs (/i,p.), (e, e) or (p, e). 
Comparing to the mode D~^ —>■ one can write 


Br 


D0^ei£2 


f3^ 


2 .2 ^^2 


mjj a 


T, 


D 


‘ D 


Ucd\^rncm^ 2 

2 + 6 ^ 


Br 


D+^/i+u 


^ 11.35 (3 


( 86 ) 


The corresponding branching ratio for the three body modes c —ulilj is given by 

0.343/32 (a2 + 62 ) / 2 . 

We have evaluated the parameters /3, a and b in several models with multiple Higgs 
scalar doublets |@,[^ and computed the branching ratios for rare decay modes of the 
D^. We find that the branching ratios for these modes can be as large as 


Broo- 




8 X 10 


-10 


Br^o. 


4 X 10 


-14 


~ 7 X 10-'" ,(87) 


with the corresponding three body modes having branching ratios smaller than these by 
about a factor of 30. While still small, these values are greatly enhanced over those in 
the SM. 


3.2.2 FCNC in Horizontal Gauge Models 

The gauge sector in the Standard Model has a large global symmetry which is broken by 
the Higgs interaction. By enlarging the Higgs sector, some subgroup of this symmetry 
can be imposed on the full SM lagrangian and break the symmetry spontaneously. This 
family symmetry can be global as well as gauged |48]. If the new gauge couplings are 


very weak or the gauge boson masses are large, the difference between a gauged or global 
symmetry is rather difficult to distinguish in practice. In general there would be FCNC 
effects from both the gauge and scalar sectors. Here we consider the gauge contributions. 

Let us construct a simple toy model as an example. Consider a family symmetry 
SU{2)h under which the left-handed quarks (where the superscripts denote the weak 
flavor eigenstates) 
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and the corresponding left-handed leptons 


/i 


transform as members oi cin Ih = 1/2 family doublet. The third family is assumed to 
have Ih = 0. The SU{2)h symmetry in this model can be thought of as a remnant of an 
SU{3)h family symmetry which has been broken to SU{ 2 ) x U{ 1 ). If {G/} are the gauge 

/ dO ' 


helds corresponding to the SU{2)h and we denote = 
then the gauge interactions are 


, i’ui = 


u 


9 




After the symmetry is broken, the mass eigenstate basis is given by 


d 

s 


= U, 




u 

c 


= u„ 


u 


e 


= Ue 


9- 


, etc, 


( 88 ) 


(89) 


The matrices Uu,Ud and each contain one angle, 6 f, and three phases. After the 
symmetry is broken, the three gauge bosons acquire different masses, m*. If the phases 
are ignored, the matrix elements for the processes of interest are: 




sin 29u cos 6 ^ 


rrin 


cos 2 du cos 2 de 

9 

mf 


cos 26u sin 26^ 


2 

mf 


■I-2 + 

mi 


/r(l + 75)/i, (90) 

sin26*„ sin26'p, „ , 

h(l +75)e 


ml 


Corresponding expressions exist for decay modes, with 9d replacing 9u- To proceed 
further, let us make the simplifying assumption that mi ~ m 2 m 3 and that the mixing 
angles are small. Then, using the constraints from the kaon system, namely the bounds 
on Kl —^ e/i and the known rate for Ki /i/2, we hnd that the branching ratios for 
charm decay modes can be as large as 


Btdo. 




3.10 


-10 


and i3r£)0^„±e=F ~ 2.10 


-13 


(91) 


which are enhanced over the SM expectations. 


3.2.3 Extra Fermions 

Additional fermions beyond those in the three families of the SM can contribute to a 
variety of rare charm decays and can serve to remove the effective GIM cancellation 
inherent to these transitions in the SM. Let us first consider the effect of an SU(2) singlet 
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down-type Q=-l/3 quark of the kind that occurs in Eq models |^. This h' quark will 
contribute in the loop diagrams which mediate decays such as For a 

mass my ~ 250 GeV, the mixing with u and c quarks given by Xy = VuyV*y is constrained 
by the h' contribution to Amo- With the current bound on xd {xd = Am^/r b) of about 
3% p, \y has to satisfy Xy < 0.003. The h' contribution to —>• can then be of 

order 

, (92) 

which is two orders of magnitude above the SM value. There will be similar enhancements 
for modes such as D —y 7iii, D —> p^^ which would be experimentally detectable. We note 
that an additional fourth family down-type quark belonging to a SU (2 )l doublet would 
have an identical effect. 

When the SM is extended by adding extra lepton doublets or extra neutral singlets, 
the decay mode —> pe can be generated (in a similar fashion as Kl ^ pe) only if there 
are non-degenerate neutrinos and nonzero neutrino mixings HI- We display the relevant 
box-diagram in Fig. [^. The associated matrix element can be written as 


Moo^iie = 2^2 uT^v 


(93) 


where B is given by |]l3 
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- E 

a,k 


afi 


U(ye ^ck^uk^a^k 


(1 
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In Xk 


x„)(l - 
In Xn 


Xk) 


Xa-Xk V (1 - 2;fc)2 (1 - X«)2 


(94) 


In the above, the greek and latin indices run respectively over the neutral leptons and 
negatively-charged quarks, Uap and Vjk are respectively mixing-matrix elements for lep¬ 
tons and quarks, and Xk = ml/M^. In the excellent approximation that ~ 0 for 


a = Ue, z/„, Ut- and x* = 0 for i = d, the expression for B becomes [52 


B — 



— InXo -t- 


InxAT 


Inxfe -t- 


(1 -X7v)2 

InxAT ^ 


(1 -X7v)2 


(95) 


for a fourth generation neutral lepton mass of ruAr ~ 50 GeV. This result varies rather 
slowly as mjv increases to larger values up to and beyond M^. The decay rate for pe 

is then given by 


^ DO ^^le — 


G]pM^fDm^,B 

47|-2 




(96) 
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Figure 15: Box diagram mediating ^ /re. 


The mixing (UNeUN/i) for mA? > 50 GeV is constrained by the limit on 
be [E^,|T6i less than 5.6 x 10and hence we infer 




to 


TdO—»/ xe 


< 8.62 X 10-27 GeV , 
< 1.3 X 1020 sec-7 _ 


(97) 


The branching ratio for —>■ /re is thus bounded by 


i3r£)0^^-e+ < 5.2 X 10 7® or i3roo^^-e++^+e- < 1.0 x 10 7"^ . (98) 

If the heavy neutral lepton iV° is an SU{2) singlet rather than a member of a doublet, 
the same result is obtained, even though the GIM suppression is absent |l52| , |5^ . Hence 
any observation of —> /re with > 10-7^ cannot be explained by mixing with 

a heavy neutrino. 


3.2.4 Extra Dimensions 


Attempts to address the hierarchy problem by exploiting the geometry of space-time have 
led to extra dimensional theories which have verihable consequences at the TeV scale. 
These theories make use of the idea that our universe lies on a (3 -|- l)-dimensional brane 
which is embedded in a higher D-dimensional space-time, D = {1 + 3 + S), known as the 
bulk. The size and geometry of the bulk, as well as the held content which is allowed to 
propagate in the bulk, varies between different scenarios. Upon compactihcation of the 
additional dimensions, all bulk helds expand into a Kaluza-Klein (KK) tower of states 
on the (3 -|- l)-brane, where the masses of the KK states correspond to the ^-dimensional 
kinetic motion of the bulk held. The direct observation or indirect ehects of the KK states 
signals the existence of extra dimensions. 


There are various potential contributions to rare decays within these scenarios: 

(i) In the case of large, hat toroidal extra dimensions |^, gravity alone propagates 
in the bulk and the resultant bulk graviton KK tower states, couple with inverse 
Planck scale strength and have very hne mass splittings given by 1/Rc ~ IQ-^eU to a 
few MeV, where Rc is the common compactihcation radius of the additional dimensions. 
They may be radiated in rare decays such as c —> u -|- and subsequently appear as 
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missing energy. The bulk graviton KK states couple to the conserved stress-energy tensor, 
giving a contribution to this process of order where Mj^ is the fundamental scale 

of gravity in the higher dimensional space and is assumed to be of order a TeV in these 
models. 


(ii) If the extra dimensions are of size TeV“^, then the Standard Model gauge helds 
may propagate in the bulk and hence expand into KK towers . The KK tower states of 
the 7 , Z, W, and gluon may participate in rare transitions in a variety of ways. However, 
precision electroweak data constrain the mass of the hrst gauge KK excitation to be in 


excess of 4 TeV |57], and hence their contributions to rare decays are small [58 


(hi) If the Standard Model fermions are localized at specihc points within a TeV“^- 
sized extra dimension, then they obtain narrow gaussian-like wave functions in the extra 
dimension with a width much smaller than the compactihcation radius. In this case, the 
fermion mass hierarchy may be explained and FCNC are suppressed by the small overlap 
of the wave functions for the different flavors. 


(iv) The last possibility is the Randall-Sundrum model of localized gravity |Q, based 
on a non-factorizable geometry in 5-dimensional Anti-de-Sitter space. In this case, the 
Standard Model gauge and matter helds, as well as gravity, are allowed to propagate in 
the warped extra dimension. The hrst bulk graviton KK excitation mass is of order a 
TeV and hence does not participate in rare decays. However, the hrst gauge and fermion 
KK excitations are lighter and may have interesting consequences in rare transitions . 
In models of this type, it is possible to generate tree-level FCNC which may produce 
observable ehects in rare charm decays. 


3.3 Strong Dynamics 

The possibility that new strong interactions are responsible for electroweak symmetry 
breaking (EWSB) and/or fermion masses has important consequences for havor physics. 
The SM with one Higgs doublet already requires the presence of new dynamics at a scale 
A in order to avoid triviality bounds. The physics above the cutoh scale gives rise to 
the scalar sector via bound states and is connected in some fashion to the the generation 
of havor. For instance, technicolor theories require extended technicolor, whereas the 
generation of the (large) top quark mass may require a top-condensation mechanism. In 
general the generation of fermion mass textures leads, in one way or another, to FCNC. 
Here we examine some of the potential ehects in rare charm decays and their relation to 
other phenomenological constraints. 
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3.3.1 Extended Technicolor 


In standard technicolor theories both fermions and techni-fermions transform under the 
new gauge interaction of Extended Technicolor (ETC). The condensation of techni-fermions 
leading to EWSB leads to fermion mass terms of the form 

m,~^{fT>ETC . (99) 

■^ETC 


The ETC interactions connect ordinary fermions with techni-fermions, as well as fermions 
and techni-fermions among themselves. The relevant sources of FCNC in technicolor 
models divide into two classes: those associated with the technicolor sector and those 
where the diagonal ETC gauge bosons acting on ordinary fermions give rise to FCNC 
through dimension-six operators. 


The first case gives rise to operators mediated by ETC gauge bosons. These in turn 
have been shown to give rise to FCNC involving the Z-boson, 


e' 


mr_ 


Snv sin 29]^ 


UcuZ^ and ^ 


2 




Stiv sin 29w 


( 100 ) 


where is the unitary matrix rotating left-handed up-type quark fields into their mass 
basis and is a model-dependent quantity of 0{1). The induced flavor-conserving Z 
coupling was hrst studied in Ref. [^] and flavor-changing effects in B decays have been 
Ml65i 


The flavor-changing vertices in Eq. 


induce contributions 


examined in Refs 

to c —*• These appear mostly as a shift in the Wilson coefficient C'io(Miy), 




nir sin^ 9 


w 


a 


~ 0.02 , 


( 101 ) 


where we make the assumption f/^ ~ A ~ 0.22 (he., one power of the Cabibbo angle) and 
we take rric = 1.4 GeV. Although this represents a very large enhancement with respect 
to the SM value of Cio(M[v), it does not translate into a large deviation in the branching 
ratio. As mentioned previously, these are dominated by the mixing of the operator O 2 
with Og, leading to a very large value of Cg^. The contribution in Eq. (p.01| ) represents 
only a few percent effect in the branching ratio with respect to the SM. On the other 
hand, the interaction in Eq. (|10(]| ) can also mediate The corresponding 

amplitude is 

- ^cu ^ ^ dwfn , (102) 

which should be compared to Eq. (^0]). This results in the branching ratio ~ 

0.6 X 10“^°, which although still small, is not only several orders of magnitude larger than 
the SM short distance contribution but also more than two orders of magnitude larger 
than the long distance estimates. 
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Finally, the FCNC vertices of the Z boson in Eq. (|100|) also give large contributions 
to c —*• uvv. The enhancement is considerable and results in the branching ratio 


^^etc 


D+—^XuUu — ^ 


uL 

0.2 


2 X 10 


-9 


(103) 


The second class of contributions from technicolor models comes from the diagonal 
ETC gauge bosons. These generate four-quark interactions which refer to a mass scale 
constrained by D^-D^ mixing to be approximately M > 100 TeV |^], thus making such 
effects very small in rare charm decays. 


3.3.2 Top-condensation Models 

Top-condensation models postulate a new gauge interaction that is strong enough to 
break the top-quark chiral symmetry and give rise to the large top mass. The various 
realizations of this basic idea have one common feature: flavor violation. Since the new 
interaction must be non-universal, it mediates FCNC at tree level. This arises because 
the mass matrix generated between the top-condensate and the other flavor physics gives 
rise to the lighter fermion masses {e.g. ETC in top color-assisted technicolor and 
is not aligned with the weak basis. Diagonalization of this mass matrix will then leave 
FCNC vertices of the so-called ‘topcolor interactions’ since they couple preferentially to 
the third generation. The exchange of top-gluons and topcolor gauge bosons will generate 
four-fermion couplings of the form 

Ana, cot 9^ (fy'^T^c) 


M2 
Arra. tan 9^ 


M2 


l/cu (h7^T“c)(c7^T“c) 


Arras 


f/e. (u^,T^c)(irTV , 


(104) 


where = (t b), Uij = and M is the mass of the exchanged color-octet gauge 

boson. The hrst term comes from rotating two top-quark helds via the strongly coupled 
topgluon, with the strong interaction being reflected in the factor cot^ 9 ~ 22. The 
second term corresponds to a topgluon which is weakly coupled to the hrst and second 
generations. In the third term, which gives the largest contribution, the topgluon couples 
strongly to the third generation quark current but weakly to the (he) current, giving rise 
to a gluon-like coupling. The one-loop insertion of the hrst and/or third terms in Eq. ( 1U4 ) 
would result in contributions to the operators Og and Oig. However, a term analogous 
to the second term in Eq. (|104|) but with the cl quark rotated to a would contribute 


to D^-D^ mixing. The current experimental bound on Amo taken from Eq. 


implies 
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that 


M 


Re[Uc, 


> 140 TeV . 


(105) 


In standard Topcolor Assisted Technicolor models, this constraint is not binding on the 
top-gluon mass since the up-sector rotation matrices are taken to be nearly diagonal |]67 


However, once it is satished, the bound of Eq. ( [105|) implies that all effects in rare charm 
decays are negligible. Similarly, this also applies to the topcolor Z' arising from the 
strongly coupled U{1)y- 


4 Conclusions 

We have extensively evaluated the potential of rare charm decays to probe physics beyond 
the SM. In Section we computed the SM rates for a variety of decay modes; incorporating 
the hrst evaluation of the QCD corrections to the short distance contributions, as well 
as a comprehensive study of long range effects. This extends our earlier work in Ref. |^, 
where we concentrated solely on radiative decays. We have shown that although, just as 
in the radiative modes, it is still true that long distance contributions dominate the rates, 
there are decay channels where it is possible to access the short distance physics. This 
is particularly true for the case oi D ^ decay modes such as D ^ and 

D —> away from the resonance contributions in the low dilepton mass region. This 

is illustrated in Figs. Q and 2, where we see that for low dilepton invariant mass the sum of 
long and short distance effects leaves a large window where physics beyond the SM can be 
observed. Although the uncertainties in our calculation of the long distance contributions 
to this mode are still sizable (roughly of 0 (1)) it is clear that at low dilepton masses 
new physics effects that are an order of magnitude or more larger than the short distance 
SM signal can be detected. This is not the case in the resonance region where the 0, a; 
and p contributions take the rates to values just below current experimental bounds, in 
a situation analogous with radiative decays such as D —> py. We compile our predictions 
for the SM rates in Table 

In Section we explored the potential of these decays to constrain new physics. In 
the case of the MSSM, we examined the sensitivity of rare charm decays to non-universal 
soft breaking in the squark mass matrices. We found that large effects are possible in 
D —^ and particularly in D ^ pi^i~, as can be seen from Figures 0 and 0. The 

effect in the vector mode is amplihed by the heightened sensitivity of this decay channel 
to the photonic penguin, which carries a large enhancement since the gluino helicity flip 
replaces the usual charm quark mass insertion. This effect, unfortunately, is obscured 
in radiative decays such as D —>■ py due to the overwhelming long range effects. It 
can therefore, only be observed by examination of the full dilepton mass spectrum in 
D —. We conclude that an important fraction of parameter space in the MSSM 
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with non-universal soft breaking can be explored if sensitivities of the order of 10 ® to 
10“^ in the kinematic region of interest are reached. 

We also considered the effects of R-parity violating conplings in supersymmetry. We 
found that the current upper limit on the decay D yields the best constraint on 

the product \ 22 k^ 2 ik Eq. ([781)). Thus rare charm decays already constrain R-parity 
violating effects! Our results are summarized in Table § for the predictions with R-parity 
violation effects, assuming the couplings saturate their current bounds. We have also 
shown that the forward-backward asymmetry for leptons Ap^ in is quite 

sensitive to these effects (c/. Figure |T^. More generally, App is negligibly small in the SM 
due to the fact that the vector coupling of leptons is enormously enhanced with respect 
to the axial-vector coupling by the presence of vector mesons. Thus, any observation of 
ApB would point to the presence of new physics. 

We also considered the effects of other non-supersymmetric extensions of the SM 
including multi-Higgs models, horizontal gauge models, a fourth generation, extra dimen¬ 
sions, as well as models with strong dynamics such as extended technicolor and topcolor. 
These scenarios give sizeable enhancements in some of the modes. 

We conclude that these rare charm decay modes are most sensitive to the effects of 
non-universal supersymmetry breaking as well as to R-parity violating couplings. It is 
then important to push for increased sensitivity of the experiments, preferably to below 
10“® in order to highly constrain these effects. This is in stark contrast with the situation 
in the radiative modes, where sensitivity below 10“® — 10“® may not illuminate short 
distance physics. The dilepton modes should be pursued by all facilities to highest possible 
sensitivity. 
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Decay Mode 

Experimental Limit 

Brs.D. 

Btl.d. 

D+ ^ X+e+e- 


2 X 10-® 


D~^ — >■ 

< 4.5 X 10-^ 


2 X lO-o 

D~^ —>• jji~ 

< 1.5 X 10"^ 


1.9 X lO-o 

D+ ^ p+e+e" 

< 1.0 X 10-^ 


4.5 X lO-o 

DO ^ X^e+e- 


0.8 X 10"® 


— >• 7r0e+e“ 

< 6.6 X 10“^ 


0.8 X lO-o 

£)0 p^e^e~ 

< 5.8 X 10-^ 


1.8 X lO-o 

DO ^ p0^+^- 

< 2.3 X 10-^ 


1.8 X lO-o 

D* ^ Xiuu 


1.2 X 10-^® 


D~^ —>■ TT^VV 



5 X 10-^0 

DO ^ 



2.4 X 10-10 

Dg —>■ TT^vV 



8 X 10-10 

DO ^ 77 


3 X 10-1^ 

few xlO-® 

DO ^ p+p- 

< 3.3 X lO-o 

10-18 

few X 10-1® 

DO ^ e+e- 

< 1.3 X 10-^ 

(2.3 -4.7) X 10-24 


DO ^ 

< 8.1 X lO-o 

0 

0 

— >• 7r’''p^e^ 

< 3.4 X 10"^ 

0 

0 

DO ^ pO^igT 

< 4.9 X 10"^ 

0 

0 


Table 6: Standard Model predictions for the branching fractions due to short and long dis¬ 
tance contributions for various rare D meson decays. Also shown are the current experimental 


limits fill , IM , pd[| ■ 
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